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Results from Charmonium Decays 
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Recent results from BES and CLEOc experiments on charmonium decays using J/tp, ip' and ip" data samples 
collected in e+e — annihilation are reviewed, including the measurement of the relative phase between strong 
and electromagnetic decays of ip' , the study of the "pir puzzle" in J/ip and ip' decays, and the search for the 
non-DZ) decays of tp". The decays of \cj produced in ip 1 radiative transition are also reviewed. These new 
results shed light on the understanding of QCD. 
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1. Introduction 

BESII 1] running at BEPC and CLEOc run- 
ning at CESR are the two detectors operating in 
the T-charm energy region. Both experiments have 
collected large data samples of charmonium decays 
including 58 M J/ip events, 14 M ip' events, and 
33 pb" 1 data around ip" peak at BESII, and 4 M 
ip' events, and 281 pb" 1 ip" events at CLEOc. To 
study the continuum background in the charmonium 
decays, special data samples at the center of mass en- 
ergy lower than the ip' mass were taken both at BESII 
(Vi = 3.650 GeV) and at CLEOc {js = 3.671 GeV), 
the luminosities are 6.4 pb" 1 and 21 pb" 1 respec- 
tively. These data samples are used for the study of 
the hadron spectroscopy, the D decay properties and 
the CKM matrix, as well as the charmonium decay 
dynamics. 

In this paper, we focus on the the extensive study 
of the "p7r puzzle" in J/ip and ip' decays, the relative 
phase between strong and electromagnetic amplitudes 
of ip' decays, and the non-Z?L> decays of ip". The study 
of Xcj decays are also reviewed. 

It should be noted that performance of the CLEOc 
detector is much better than the BESII detector, espe- 
cially in the photon detection, this makes the 4 M ip' 
events data sample at CLEOc produces results with 
similar precision as from 14 M ip' events from BESII. 



2. Relative Phase in r' Decays 

It has been determined that for many two-body 
exclusive J/ip decays, like vector pseudoscalar (VP), 
vector vector (VV), and pseudoscalar pseudoscalar 
(PP) meson decays and nucleon anti-nucleon (NN) 
decays, the relative phases between the three-gluon 
and the one-photon annihilation amplitudes are near 
90° 0, S H IE II 11 For V' decays, the available 
information about the phase is much more limited be- 
cause there are fewer experimental measurements. It 
has been argued that the relative phases xnxjp 1 decays 
should be similar to those in J / ip decays |3j, |9j , but the 
analysis of ip' to VP decays in Ref. 3] indicates this 
phase is likely to be around 180°. Another analysis of 



this mode though shows the relative phase observed 
in J/ip decays could also fit these decays ^jj, but it 
could not rule out the 180° possibility due to the big 
uncertainties in the experimental data. 

BES measured the branching fraction of ip' — > 
"13 to be (5.24 ± 0.47 ± 0.48) x 10" 5 , to- 
gether with the known branching fractions of ip' — > 
7r + 7r" and ip' — > K + K~, two possible solutions of 
the phase are found, which is either —(82 ± 29)° or 
+(121 ± 27)°, following the procedure developed in 
Ref. 12\. Benefits from the good detector perfor- 
mance, CLEO measures all three pseudoscalar me- 
son pair decay modes 0. The branching fraction 
of ip' -> K° S K° L is (5.8 ± 0.8 ± 0.4) x 10" 5 , that of 
ip' -> K+K- is (6.3 ± 0.6 ± 0.3) x 10" 5 ; the signal of 
ip' — > 7r + 7T" is not significant, and the upper limit of 
the branching fraction is determined to be 2.1 x 10 -5 
at 90% C.L. CLEO measures the relative phase to be 
(95 ± 15)°, assuming there is no interference between 
continuum amplitude and the resonance decay ampli- 
tudes, which is questionable since the two components 
of the resonance decays has a relative phase around 
90°, as measured by both BESII and CLEOc exper- 
iments. Following the same procedure developed in 
Ref. fl2^. CLEO measurements also result in two so- 
lutions for the phase, one at around —80 degrees, while 
another at around +120 degrees, in good agreement 
with the BES result. The error is large due to the 
small statistics of the data samples, and there is no 
way to determine which solution is the physical one. 



3. "p7r Puzzle" in J/ip> and ip' Decay 

From perturbative QCD (pQCD), it is expected 
that both J / ip and ip' decaying into light hadrons are 
dominated by the annihilation of cc into three glu- 
ons or one virtual photon, with a width proportional 
to the square of the wave function at the origin |l4|. 
This yields the pQCD "12% rule", 
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A large violation of this rule was first observed in de- 
cays to pn and K* + K~ + cc. by Mark II ^5(, known 
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as the pn puzzle, where only the upper limits on the 
branching fractions were reported in ip' decays. Since 
then, many two-body decay modes of the ip' have been 
measured by the BES collaboration and recently by 
the CLEO collaboration; some decays obey the rule 
while others violate it [II E E El d • 



3.1. V — ► pn' Current Status 

The pn mode of the vector charmonia decays is es- 
sential for this study, since this is the first puzzling 
channel found in J/ip and ip' decays. The new mea- 
surements, together with the old information, show us 
a new picture of the charmonium decay dynamics. 



(2.4t°j±0.2) x 1(T 5 , the difference is large. Although 
a big difference exists between the BESII and CLEOc 
results, it does mean that the ip' — > pn signal exists, 
rather than the signal is completely missing as has 
been guessed before. If we take a weighted average ne- 
glecting the difference between the two measurements, 
we get B(ip' -► pn) = (3.1 ± 0.7) x 1CT 5 . 

Comparing B(ip' — > pn) with B(J/ip — > pit), one 
gets 

The suppression compared to the 12% rule is obvious. 



3.1.1. J/tp -> TT+TT^TT 

BESII measured the J/ip — > n + n~n° branching 
fraction using its J/ip and ip 1 — > J/ipn + n~ data sam- 
ples E3> an( i BABAR measured the same branching 
fraction using J/ip events produced by initial state 
radiative (ISR) events at s/s = 10.58 GeV ;2fJ|. To- 
gether with the measurement from Mark-II 15], we 
get a weighted average of B(J/ip — > n + n~n°) = 
(2.00 ±0.09)%. 

To extract the J/tp — ► pn branching fraction, 
partial wave analysis (PWA) is needed to consider 
the possible contribution from the excited p states, 
the only available information on the fraction of pn 
in J/ip — > n + n~n° was got from Mark-Ill. Us- 
ing the information given in Ref. J2l|, we estimate 

= L17 ( 1 ± 10% )' with the error from 

an educated guess based on the information in the 
paper since we have no access to the covariant matrix 
from the fit showed in the paper. From this num- 
ber and B(J/ip — > 7r + 7r~7r°) got above, we estimate 
B{J/ip -► pn) = (2.34 ± 0.26)%. This is substan- 
tially larger than the world average listed by PDG , 
which is (1.27 ±0.09)%. 

3.1.2. ip' -> n+n-n° 

ip' — > pn was studied both at BESII ^3 and 
CLEOc [lSj . After selecting two charged pions and 
two photons, clear n° signals are observed in the two 
photon invariant mass spectra, the numbers of signals 
are found to be 229 and 196 from BESII and CLEOc 
respectively, and the branching fraction of ip' — > 
n+n-n° is measured to be (18.1 ± 1.8 ± 1.9) x 10" 5 
and (18.8±i;| ± 1.9) x 10" 5 at BESII and CLEOc re- 
spectively. The two experiments give results in good 
agreement with each other. 

To extract the branching fraction of ip' — > pn, how- 
ever, BESII uses a PWA including the high mass 
p states and the interference between them, while 
CLEOc counts the number of events by applying a 
p mass cut. The branching fraction from BESII is 
(5.1 ± 0.7 ± 1.1) x 10" 5 , while that from CLEOc is 



3.2. Other Studies 

There are many more new measurements on ip' de- 
cays for the extensive study of the 12% rule (l^lialill 
HsI l23l ] , among which the VP modes are measured as 
first priority. The ratios of the branching fractions are 
suppressed for almost all these VP modes compared 
with the 12% rule; while the PP modes, K° S K° L and 
K + K~ , are enhanced compared with the 12% rule. 
The multi-hadron modes and the baryon-antibaryon 
modes are either suppressed, or enhanced, or normal, 
which are very hard to be categorized simply. The var- 
ious models, developed for interpreting specific mode 
can hardly find solution for all these newly observed 
modes. 

One observation is that many of the attempts to in- 
terpret the pn puzzle are based on the potential mod- 
els for the charmonium which were developed more 
than 20 years ago, as the B-factories discovered many 
new charmonium states |24| which are hard to be ex- 
plained in the potential models, it may indicate even 
the properties of J/ip and ip' are not as expected from 
the potential models. The further understanding of 
the other high mass charmonium states may shed light 
on the understanding of the low lying ones. 



4. Uon-DD Decays of ip' 1 



Since ip" is above the DD threshold, it decays pre- 
dominantly into charmed mesons, however, since the 
old measurements may indicate big contribution of ip" 
charmless decays [2^ , and large fraction of charmless 
decays of ip" is expected if ip" is a mixture of S- and D- 
wave charmonium states and the mixing is responsible 
for the 12% rule violation in J/ip and ip' decays |25| . 
both BESII and CLEOc experiments tried to search 
for the non-DZ) decays of ip" in both exclusive modes 
and inclusive measurement. 
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4.1. Exclusive Decays 

It has been pointed out that the continuum am- 
plitude plays an important role in measuring ip" de- 
cays into light hadrons [2{|. In fact, there are two 
issues which need to be clarified in ip" decays, that 
is whether ip" decays into light hadrons really exist, 
and if it exists, how large is it. The searches for the 
uon-DD decays are performed by comparing the cross 
sections on and off the ip" resonance peak. 

4.1.1. ip" -» tt+ty-tt 

By comparing the cross sections of e + e~ — > n + ir~n° 
at the ip" resonance peak (yG = 3.773 GeV) and at 
a continuum energy point (-^s = 3.650 GeV at BESII 
and 3.671 GeV at CLEOc) below the ip' peak, both 
BESII and CLEOc found that a(e + e~ -> tt+tt-tt ) at 
continuum is larger than that at ip" resonance peak. 
The average of the two experiments [2(| [27I are 

a{e + e- n + n-n a ) on = 7.5 ± 1.2 pb, 

er(e+e~ -> 7r + 7r~7r°) off = 13.7 ± 2.6 pb. 

The difference, after considering the form factor varia- 
tion between 3.650 and 3.773 GeV, is still significant, 
and it indicates that there is an amplitude from ip" 
decays which interferes destructively with the contin- 
uum amplitude, and makes the cross section at the ip" 
peak smaller than the pure contribution of continuum 
process. 

For the pn mode, BESII can only give upper limit 
of its cross section due to the limited statistics of 
the data sample, th e up per limit at 90% C. L. is 
found to be 6.0 pb |2|j at the ip" peak, which is 
in consistent with the measurement from CLEOc us- 
ing a much larger data sample: cr(e + e~ — > p7r)on = 
4.4 ± 0.6 pb 271; while the cross section at the con- 
tinuum is 8.0ti;4 ± 0.9 pb measured by CLEOc. 

To extract the information on the ip" — > pir branch- 
ing fraction, BESII developed a method based on the 
measured cross sections at ip" resonance peak and at 
the continuum [2(|. By neglecting the electromag- 
netic decay amplitude of ip", there are two ampli- 
tudes which contribute to the cross section at the ip" 
peak, the strong decay amplitude of ip" and the con- 
tinuum amplitude. Taking the continuum amplitude 
as a real number, the ip" strong decay amplitude is de- 
scribed as one real number for the magnitude, and one 
phase between ip" strong and electromagnetic decays 
to describe the relative phase between the two ampli- 
tudes. Since only two measurements are available (at 
ip" peak and at continuum), one can only extract ip" 
decay branching fraction as a function of the relative 
phase. BESII measurement on the upper limit of the 
e + e~ — > pw cross section at ip" peak, together with the 
CLEOc measurement of the continuum cross section 
restrict the physics region of the branching fraction 



and the relative phase as shown in Fig. ^ From the 
Figure, we see that the branching fraction of ip" — » pn 
is restricted within 6 x 10~ 6 and 2.4 x 10~ 3 , and the 
phase is between -150° and -20°, at 90% C.L. 




-180 -160 -140 -120 -100 -80 -60 -40 -20 
Phase between Vf" Strong and EM decays 

Figure 1: Physics region on B(ip" — > pn) and the relative 
phase (0) between ip" strong and electromagnetic decays 
from BESII. 

The observation of the e + e~ — > pir signal at ip" 
peak and the measurement of the cross section 27] 
at CLEOc further make the physical region in the 
branching ratio and relative phase plane smaller: the 
CLEOc measurement gives a similar outer bound of 
the physical region as BES gives, while also indicates 
the central part of the physical region in Fig. is not 
allowed by physics. By using a toy Monte Carlo to 
simulate the CLEOc selection criteria and the inter- 
ference between the resonance and continuum ampli- 
tudes, we found that the ip" — > pir branching fraction 
could be either (2.1 ± 0.3) x 10" 3 or (2.4t 3 ;*) x 10~ 5 
from the CLEOc measurements, if the relative phase 
between ip" strong and electromagnetic decay ampli- 
tudes is —90° as observed in J/ip and ip' decays [lCj . 

4.1 .2. Other Exclusive Modes 

The ip" decays into light hadrons were searched for 
in various ip" decay modes, including two-bady and 
multi-hadron modes [27], Hg, |2j|. However, only the 
comparison between the cross sections at continuum 
and those at ip" resonance peak are given, instead of 
giving the ip" decay branching fractions. In current 
circumstances, it is still not clear whether the ip" de- 
cays into light hadrons with large branching fractions, 
since, as has been shown in the pn case, there could 
be two solutions for the branching fraction, and the 
two values could be very different. 

As the luminosity at the ip" peak is large enough, 
current study is limited by the low statistics at the 
continuum: at CLEOc, the luminosity at continuum is 
more than an order of magnitude smaller than that at 
peak, this prevents from a high precision comparison 
between the cross sections at the two energy points. 
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One conclusion we can draw from the existing data 
is that the measurements do not contradict with the 
assumption that the relative phase between ip" strong 
and electromagnetic decay amplitudes is around —90°, 
and the if)" decays into light hadrons could be large. 



4.2. Inclusive Measurements 

The inclusive non-DD decays of ip" is searched 
for by measuring the DD cross section and the to- 
tal hadronic cross section above the uds continuum 
contribution at the if)" resonance peak. 

BES and CLEO measure the DD cross section |3CJ, 
l3ll | using both single tag and double tag methods, the 
results are shown in Fig. [3 Good agreement between 
BES and CLEO on the D+D~ and D°D°, as well as 
the total DD cross sections are found. The weighted 
average of the two experiments is (6.32+Jj 'U) nb for 
the total DD cross section. 
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Figure 2: Measurements of the DD cross section from 
BES and CLEOc, and the weighted average of the 
measurements. 

The total hadronic cross section of if>" decay is 
obtained by comparing the cross sections at the ip" 
peak and the data points below the resonance peak. 
The contributions of the radiative tails of J I if) and 
ip' are also subtracted. By comparing this cross sec- 
tion and the DD cross section, BES found a non-DD 
decay branching fraction of (14.0 ±_1.7 ± 6.0)% |H, 
while CLEO measured the non-DD cross section of 
(-0.01 ± 0.08±g;3£) nb [H, which corresponds to an 
upper limit of the non-DD decay branching fraction 
of 11% at 90% C. L. 

Although the BES and CLEO measurements are 
not in contradiction considering the large uncertain- 
ties, the big difference between the two central values 
calls for a better understanding of the procedure from 
which the total inclusive cross section is obtained. An 
obvious difference between the two measurements is 
how to consider the interference between the reso- 
nance and the continuum amplitudes. 



5. Xcj Decays 

Since each XcJ (J = 0, 1, 2) is produced about 10% 
of the ip' decays, they are studied at BES to under- 
stand the P-wave charmonium decay dynamics, as 
well as the light hadron spectroscopy. 

5.1. xcj ~^VV 

Two modes, uiui [3j| and <f)(f>, were measured re- 
cently, while the former is the first observation, the 
latter improves the precision. The results are sum- 
marized in Table [IJ together with the measurement 
of XcJ -> if* (892)^* (892) 0. These results are 
used to predict the decay branching fractions of \ c j 
to other vector meson pairs, like pp and ivcp j^, large 
double OZI suppressed amplitude is expected. 



Table I Branching fractions of \ c j — > VV (in 10 3 ), the 
results on <fxf> are BES preliminary. 



Decay mode 




Xc2 


LOLO 


2.29 ±0.58 ±0.41 


1.77 ±0.47 ±0.36 


# 


0.94 ±0.21 ±0.14 


1.48 ±0.26 ±0.23 


K* (892)iT (892) 


1.78 ±0.34 ±0.34 


4.86 ± 0.56 ± 0.88 



5.2. xcj - 

The importance of the Color Octet Mechanism 
(COM) for Xcj decays has been pointed out for many 
years |36j , and theoretical predictions of two-body ex- 
clusive decays have been made based on it. Recently, 
new experimental results on Xcj exclusive decays have 
been reported pV?! l38l | . COM predictions for many Xcj 
decays into meson pairs are in agreement with exper- 
imental values, while predictions for some decays into 
baryon pairs, for example, the branching fractions of 
Xcj — * A A, disagree with measured values. For fur- 
ther testing of the COM in the decays of the P-wave 
charmonia, measurements of other baryon pair decays 
of XcJ, such as XcJ —* is desired. 

The measurement of Xco ~ * S _ S + is helpful for un- 
derstanding the Helicity Selection Rule (HSR) [sflj ]. 
which prohibits Xco decays into baryon antibaryon 
{BB) pairs. However, the measured branching ratios 
for Xco decays into pp and AA do not vanish, demon- 
strating a strong violation of HSR in charmonium de- 
cays. Measurements of Xco decays into other baryon 
anti-baryon pairs would provide additional tests of the 
HSR. 

The measured branching fractions or upper limits 
are summarized in Table [H] [4Tl |. along with some 
theoretical predictions. Theoretically, the quark cre- 
ation model (QCM) predicts B(xcO — > = 
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(2.3 ± 0.7) x 1CP 4 , which is consistent with the ex- 
perimental value within la. For Xci & n d Xc2 decays 
into S _ S + , the measured upper limits cover both the 
COM and QCM predictions. Within 1.8cr the branch- 
ing fraction of Xco ~ * does not vanish. For fur- 
ther testing of the violation of the HSR in this decay, 
higher accuracy measurements are required. 

Table II The comparison of the branching fractions or 
upper limits for \ c j — > E~H + between experimental 
values and theoretical predictions. The COM predictions 
are from Ref. |4l|. and the quark creation model (QCM) 
predictions are from Ref. |42|. The numbers are in unit 
of 10 -4 , the upper limits are at 90% C. L. 



Channel 


Branching Fraction 


COM 


QCM 




5.3 ±2.7 ±0.9 




2.3 ±0.7 




or < 10.3 






Xcl 


< 3.4 


0.24 




Xc2 


< 3.7 


0.34 


0.48 ±0.21 



5.3. xcj - PPP 

Xco decays into three pseudoscalars is forbidden by 
the spin-parity conservation, while \ci decays is sup- 
pressed due to high orbital angular momentum. BES 
tries to measure the branching fractions of Xci decays 
into KgK + ir~ ± c.c. and i]ir + ir~. Significant signals 
are observed and the branching fractions are measured 
as 

B(xci -> K%K + n~ + c.c.) = (4.1 ± 0.3 ± 0.7) x 10 -3 , 

B(Xd -> = (6.1 ± 0.8 ± 1.0) x 10" 3 . 

The KgK + ir~ + c.c. events are mainly produced via 
if* (892) intermediate state, and the nir + Tr~~ events via 
/ 2 (1270)J7 or ao(980)7r. Except xa -> K%K+tt- ± 
c.c, all other modes are first observations. Xc2 - > 
KgK + ir~ +c.c. is also observed for the first time with 
a branching fraction of (0.8 ± 0.3 ± 0.2) x 10~ 3 . These 
results from BESII experiment are preliminary. 

5-4. XcO -> SS 

Partial wave analysis of Xco ~ * n + n~ K + K~ is per- 
formed [43] using Xco produced in -0' decays at BESII 
to study the pair production of scalars. In 14 M pro- 
duced if)' events, 1371 ip' — > jXcO, XcO — * Tr + n~K + K~ 
candidates are selected with around 3% background 
contamination. 

Besides (tttt)(KK) and (Kir)(Kir) modes which are 
used to study the scalars, (Kmr)K mode which leads 
to a measurement of #i (1270).^ and K X (\4QQ)K de- 
cay processes is also included in the fit. The PWA 



results are summarized in Table IIIII From these re- 
sults, we notice that scalar resonances have larger de- 
cay fractions compared to those of tensors, and such 
decays provide a relatively clean laboratory to study 
the properties of scalars, such as /o(980), /o(1370), 
/o(1710), and so forth. The upper limits of the pair 
production of the scalar mesons which are less signif- 
icant are determined at the 90% C.L. to be at 10 -4 
level. 



Table III Summary of the XcO — > tv + tt~ K + K~ results, 
where X represents the intermediate decay modes, and 
s.s. indicates signal significance. 



Decay mode 
(X) 


B[XcO -» X -» 
-K+TT-K+K-] (10 -4 ) 


s.s. 


/o(980)/ (980) 


3.46 ± 1.08^1.5? 


5.3a 


fo (980) f (2200) 


8.42 ± 1.42+2 ^ 


7.1a 


/o(1370)/ (1710) 


7.12 ± 1.85±?:|| 


6.5a 


K* (892)° K* (892)° 


8.09 ± 1.69l?;gg 


7.1a 


A" * (1430)7?,* (1430) 


10.44 ± 2.37l?;go 


7.2a 


^(1430)^1(1430) + c.c. 


8.49 ± 1.66±i |g 


8.7a 


ifi(1270) + 7f- ±c.c, 






#1(1270) -> Kp(770) 


9.32 ± 1.83ll;g4 


8.6a 


Ki(U00) + K~ +c.c, 






if i (1400) -> if*(892)7r 


< 11.9 (90% C.L.) 


2.7a 



The above results supply important information 
on the understanding of the natures of the scalar 
states |35l , as well as the decay dynamics of Xcj decays 
into pair of scalar particles. 



6. Summary 

There are many new results in charmonium de- 
cays from BESII and CLEOc experiments. The decay 
properties of the vector charmonium states have been 
studied for more than three decades, but they are still 
far from being understood, one extreme example is the 
"pTr puzzle" in J/tjj and ip' decays. Further studies of 
all these are expected from the BESIII at BEPCII 
which will start its data taking in late 2007. 
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